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ABSTRACT 

Aims. We investigate wave propagation and energy transport in magnetic elements, which are representatives of small scale magnetic 
flux concentrations in the magnetic network on the Sun. This is a continuation of earlier work by Hasan et al. (2005). The new 
features in the present investigation include a quantitative evaluation of the energy transport in the various modes and for dilferent 
field strengths, as well as the effect of the boundary-layer thickness on wave propagation. 

Methods. We carry out 2-D MHD numerical simulations of magnetic flux concentrations for strong and moderate magnetic fields for 
which y6 (the ratio of gas to magnetic pressure) on the tube axis at the photospheric base is 0.4 and 1.7, respectively. Waves are excited 
in the tube and ambient medium by a transverse impulsive motion of the lower boundary. 

Results. The nature of the modes excited depends on the value of p. Mode conversion occurs in the moderate field case when the 
fast mode crosses the y6 = 1 contour. In the strong field case the fast mode undergoes conversion from predominantly magnetic to 
predominantly acoustic when waves are leaking from the interior of the flux concentration to the ambient medium. We also estimate 
the energy fluxes in the acoustic and magnetic modes and find that in the strong field case, the vertically directed acoustic wave fluxes 
reach spatially averaged, temporal maximum values of a few times 10* erg cm"^ s"' at chromospheric height levels. 
Conclusions. The main conclusions of our work are twofold: firstly, for transverse, impulsive excitation, flux tubes/sheets with strong 
fields are more efficient than those with weak fields in providing acoustic flux to the chromosphere. However, there is insufficient 
energy in the acoustic flux to balance the chromospheric radiative losses in the network, even for the strong field case. Secondly, the 
acoustic emission from the interface between the flux concentration and the ambient medium decreases with the width of the boundary 
layer. 

Key words. Sun: magnetic field - Sun: photosphere ~ Sun: faculae, plages - Magnetohydrodynamic (MHD) - Waves 
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1. Introduction 

Quantitative studies of wave propagation in magnetically struc- 
tured and gravitationally stratified atmospheres help to identify 
various physical mechanisms that contribute to the dynamics of 
the magnetic network on the Sun, and to develop diagnostic tools 
for the helioseismic exploration of such atmospheres. Magnetic 
fields play an important role in the generation and propagation of 
waves. The aim of this work is to attempt a better understanding 
of this process in the magnetized solar atmosphere. We have car- 
ried out a number of numerical simulations of wave propagation 
in a two-dimensional gravitationally stratified atmosphere con- 
sisting of individual magnetic flux concentrations representative 
of solar magnetic network elements of different field strengths. 

While the magnetic field in the internetwork regions of the 
quiet Sun is mainly shaped by the convective-granular flow with 
a predominance of horizontal fields and rare occurrence of flux 
concentrations surpassing 1 kG, the magnetic network shows 
plenty of flux concentrations at or surpassing this limit with a 
typical horizontal size-scale in the low photosphere of 100 km. 
These "magnetic elements" appear as bright points in G-band 
images near disk center and they can be well modeled as mag- 
netic flux tubes and flux sheets. Their magnetic field is mainly 
vertically directed and they are in a highly dynamical state 
( |Muller|1985]|Muller et al.[T994l|Berger & Title|1996|[200T] i. 



Different from the shock induced Can H2v and bright 
points in the cell interior, the network in the chromosphere is 
seen to be continuously bright ( [Lites et al.||1993 



Sheminova 



et al.||200 5), which asks for a steady heating mechanism. It is 
also seen that the Ca ii H and K line profiles from the network 
are nearly symmetric ( Grossmann-Doerth et al.|1974| l. 



Several numerical investigations have been carried out to ex- 
plain these observations. Early works modelled the network as 
thin flux tubes and studied the transverse and longitudinal waves, 
which can be supported by them, excited by the impact of gran- 
ules. These works failed to explain the persistent emission that 
was seen in observations of the Can H and K lines. When high 
frequency waves, generated by turbulence in the medium sur- 
rounding flux tubes, were added (Hasan et al. 2000 1, the observa- 



tional signature of the modelled process became less intermittent 
and was in better agreement with the more steady observed in- 
tensity from the magnetic network. Later works examined mode 
coupling between transverse and longitudinal modes in the mag- 
netic network, using the nonlinear equations for a thin flux tube 
(eg. [Hasan & Ulmschneider 2004j). All these studies modelled 
the network as consisting of thin-flux tubes, an approximation 
that becomes invalid at about the height of formation of the emis- 
sion peaks in the cores of the H and K lines. Also, this approx- 
imation does not treat the dispersion of magnetic waves caused 
by the variation of the magnetic field strength across the flux 
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concentration and it does not take into account the emission of 
acoustic waves into the ambient medium. 



Numerical simulations by [Rosenthal et al. (20021 and 



Bogdan et al. (2003J, studied wave propagation in two- 



dimensional stratified atmospheres in the presence of a magnetic 
field. They recognized and highlighted the role of refraction of 
fast magnetic waves and the role of the surface of equal Alfven 
and sound speed as a wave conversion zone, which they termed 



the magnetic canopy. While the thick flux sheets of Rosenthal 
et al. (20021 and Bogdan et al. ( 2003 [ l were a more realistic 



model for the network, they also assumed that the magnetic field 
was potential. Considering that the gas pressure, kinetic energy 
density, and the energy density of the magnetic field are all of 
similar magnitude in the photosphere, this assumption is proba- 
bly not satisfied. 

Cranmer & van Ballegooijen ( 2005| l modelled the network 



as consisting of a collection of smaller flux tubes that are spa- 



tially separated from one another in the photosphere. Hasan et al. 
( |2005 ) performed MHD simulations of wave generation and 
propagation in an individual magnetic flux sheet of such a col- 
lection and confirmed the existence of magneto-acoustic waves 
in flux sheets as a result of the interaction of these magnetic flux 
concentrations with the surrounding plasma. They used a non- 
potential field to model the network. They speculated that a well 
defined interface between the flux sheet and the ambient medium 
may act as an efficient source of acoustic waves to the surround- 
ing plasma. In a later paper, jHasan & van Ballegooijen ( 2008| l 
showed that the short period waves that are produced as a result 
of turbulent motions can be responsible for the heating of the 
network elements. 

[Cally (2005, 2007) provided magneto-acoustic-gravity dis- 
persion relations for waves in a stratified atmosphere with a ho- 
mogeneous, inclined magnetic field and discussed the process 
of mode transmission and mode conversion. Khomenko et al.l 
( |2008| l presented results of nonlinear, two-dimensional, numer- 
ical simulations of magneto-acoustic wave propagation in the 
photosphere and chromosphere in small-scale flux tubes with 
internal structure. Their focus was on long period waves with 



periods of three to five minutes. 'Steiner et al. 



( 2007]) considered 



magnetoacoustic wave propagation in a complex, magnetically 
structured, non-stationary atmosphere. They showed that wave 
travel-times can be used to map the topography of the surface 
of thermal and magnetic equipartition (/? - 1) of such an at- 
mosphere. Hansteen et al.| ^2006 ) and |De Pontieu et al. ( |2007| l 
performed two-dimensional simulations covering the solar at- 
mosphere from the convection zone to the lower corona. They 
showed how MHD waves generated by convective flows and os- 
cillations in the photosphere turn into shocks higher up and pro- 
duce spicules. 

Despite these efforts, the physical processes that contribute 
to the enhanced network emission are still not well understood. It 
is well known, that small scale magnetic elements have varying 



field strengths, ranging from hecto-gauss to kilo-gauss ( |Solanki 
[1993] [Berger et aLl[2^4l ). This suggests that the yS = 1 layer 
in these elements varies considerably in height, which in turn 
should affect the wave propagation in them ( Schaff'enberger e t aT] 
[2005 ). Hasan et al. (2005 ) and H asan & van Ballegooijen (200^ 
argue that the network is heated by the dissipation of magnetoa- 
coustic waves. However, these works did not provide quantita- 
tive estimates of the energy flux caiTied by the waves. This is 
the main focus of the present investigation, where we examine 
wave propagation in magnetic elements with different magnetic 
field strengths. We also study the effects of varying the interface 



thickness between the flux sheet and the ambient medium on the 
acoustic wave emission in the ambient medium. 

The outline of the paper is as follows. Section |2j discusses 
the construction of the initial equilibrium model and Sect.[3]the 
boundary conditions and method of solution for the simulation. 
In Sect. |4] the dynamics and in Sect. |5] the energetics is dis- 
cussed. Section[6] discusses the effects of boundary layer thick- 
ness on the acoustic wave emission. Sect. |7] summarizes the re- 
sults and Sect. |8] contains the main conclusion and a discussion 
of the results. 



2. Initial equilibrium model 

The initial atmosphere containing the flux sheet is computed in 
cartesian coordinates using the numerical methods described in 
[Steiner et a"Ll ( |1986 ) (see also [Steiner et al.,,2007 ). The method 
consists of initially specifying a magnetic field configuration and 
the pressure distribution in the physical domain. The magnetic 
field can be written in terms of the flux function tfr(x, z) as 



B, 



diff 



dx 



(1) 



The gas pressure as a function of height and field line (flux 
value), p{^, z), is given by. 



Pp+Pc 

Po 

Pp+Po 



(.Po + P2>P') if < lA < lAi. 
(fl((A-^i)"+M'A-<Ai)^ 



(2) 

+c(i// - ij/i) + d) if ifr[<i^<ip2. 



Pp 



Pc, 



if 1^2 < tA < lAn 



where the constants a, b, c, and d are chosen such that the pres- 
sure and its first derivative with respect to (/r is a continuous func- 
tion of i// and where we choose n = 8. The three cases in Eq. (j2]l 
refer to the inteiior of the flux sheet, its boundary layer, and the 
ambient medium from top to bottom, respectively. The quantities 
Pp, pc and pq are defined as. 



Pp = Po.pe 



Pc = Po,ce 

Po - PO,p + PO.c, 



(3) 
(4) 
(5) 



where and are the photospheric and chromospheric pres- 
sure scale heights, respectively. We choose Hj, = 110 km and 
He = 220 km. Bq is the magnetic field strength on the axis of the 
flux sheet at the reference height z = 0. The equation of motion 
along each field line reduces to the hydrostatic equation 



B-[Vp-pg]^0, 



(6) 



which defines us the density distribution and with it the temper- 
ature field. From the equation of motion perpendicular to B we 
obtain the electric current density 



7= -^Bx[V/?-pg], 

which after some manipulation reduces to 
dp 



Jy = 



(7) 



(8) 
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Physical quantity 


Moderate field 


Strong field 


Sheet Axis 


Ambient medium 


Sheet Axis 


Ambient medium 


Temperature [K] 


9142 
4758 


9142 
4758 


9142 
4758 


9142 
4758 


Density [g cm"^^] 


1.2 X 10-'2 
1.4 X 10-^ 


2.4 X 10^'- 
2.7 X 10^' 


1.2 X 10^'^ 
1.4 X 10-'' 


6.0 X 10"'^ 
6.7 X lO^'' 


Pressure [dyn cm"-] 


1.3 

4.2 X 10-* 


2.7 
8.3 X 10* 


1.3 
4.2 X 10* 


6.6 
20.6 X 10" 


Sound speed [km s"'] 


13.5 
7.1 


13.5 
7.1 


13.5 
7.1 


13.5 
7.1 


Alfven speed [km s"'] 


304 
6 


212 
0.9 


582 
12 


259 
0.6 


Magnetic field [G] 


119 
801 


117 
17 


229 
1601 


225 
18 


Plasma j3 [-] 


2.0 X 
1.7 


5.0 X 10-^^ 
7.4 X 10^^ 


6.0 X lO^-* 
0.4 


3.0 X 10-3 
1.6 X lO"* 



2000 



1S00 



1 

<3 



soo 




Fig. 1. Vertical component of the magnetic field at the base of 
the flux sheet, z=0. Red and blue curves correspond to field con- 
figurations with a sharp and a wide interface to the weak-field 
surroundings, respectively. Each configuration is subdivided into 
a case of moderate field-strength with max(B,) - 800 G and a 
case of strong field with max(B,) - 1600 G. 



z = 0) of 800 G and 1600 G, on the axis of the sheet. For each of 
these cases we consider two boundary widths: 20 km and 80 km 
at the reference height z = 0. These widths can be obtained by 
choosing appropriate values of ^ii and (/'2 in Eq. (j2]i. 

The vertical component of the magnetic field at z = is 
shown in Fig.[T]for the strong and moderate field cases. For the 
sharp interface (red curves) the vertical component of the mag- 
netic field component drops sharply, whereas in the case of a 
wide interface (blue curves) the field decreases smoothly. 

The characteristic properties of the two models are summa- 
rized in Table [T] The numbers in the first row of each quan- 
tity corresponds to the top boundary (z - 1600 km) and the 
numbers in the second row corresponds to the bottom bound- 
ary (z = km). The temperature increases monotonically from 
4758 K in the photosphere to 9142 K in the chromosphere cor- 
responding to the sound speed variation from 7.1 kms ' to 
13.5 kms '. The density and pressure at the axis of the sheet 
is the same for both the cases. We should mention that the am- 
bient magnetic field is weak (of the order of few tens of Gauss). 
As we go higher up in the atmosphere the flux sheet expands 
and becomes uniform near the top with a average field strength 
of 1 18 G and 227 G for the moderate and strong field cases, re- 
spectively. The plasma-yS on the sheet axis is 1 .69 and 0.42 at the 
base for the moderate and strong field cases. 



The new magnetic field configuration can be calculated from the 
current density using the Grad-Shafranov equation. 



5V 



(9) 



The above elliptic partial differential equation can be solved 
using standard numerical methods with appropriate boundary 
conditions. In practice we solve Eq. (|9]l on a computational do- 
main that consists of only half of the flux sheet of horizontal and 
vertical extensions of 640 km and 1600 km, respectively. This 
domain is discretized on a equidistant rectangular mesh with 
a spacing of 5 km. The left side of the domain corresponds to 
the axis of the flux sheet. The value of is prescribed at the 
left and the right side boundaries. At the top and bottom bound- 
aries we use the Neumann condition difr/dz - 0, assuming that 
the horizontal field component vanishes at these two boundaries. 
We treat two different cases corresponding to field strengths (at 



3. Method and boundary conditions 

Waves are excited in the equilibrium magnetic field configura- 
tion through a transverse motion of the lower boundary (similar 
to Hasan et al. 2005). The system of MHD equations, given in 
conservation-law form for an inviscid adiabatic fluid, is solved 

T994|l. These 



according to the method described in Steiner et al. 



are the continuity, momentum, entropy, and the magnetic induc- 
tion equations. The unknown variables are the density, p, the mo- 
menta, pV V and pV,, where Vj and V. are the horizontal and verti- 
cal components of the velocity, the entropy per unit mass, s, and 
the magnetic field, B ^ and B.. The equation of state is that for the 
solar mixture with a constant mean molecular weight of 1.297. 
For the numerical integration, the system of MHD equations are 
transformed into a system of discrete finite volume equations. 
The numerical fluxes are computed based on the flux-corrected 
transport (FCT) scheme of |Oran & Boris|fl987| i. For the indue- 
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tion equation we use a constrained transport scheme ( |Devore| 
1991| l, which automatically keeps V B = 0. The time integration 
is explicit. The scheme is of second order accuracy. 

Transmitting conditions apply to the side boundaries set by 
constant extrapolation of the variables from the physical domain 
to the boundary cells. Constant extrapolation also applies to the 
horizontal component of the momentum at the top and bottom 
boundary and to the vertical component at the bottom boundary. 
The density in the top boundary cells is determined using linear 
log extrapolation, while at the bottom boundary hydrostatic ex- 
trapolation applies. For the temperature constant extrapolation is 
used. The horizontal component of the magnetic field at the top 
and bottom boundaries are set equal to the corresponding values 
at the preceding interior point. The vertical component of the 
magnetic field is determined by the condition V ■ B = 0. 

The transverse velocity at z = is specified as follows: 



Vq sin(27rf/P) for < r < f /2 , 
for > f > P/2 , 



(10) 



where Vq denotes the amplitude of the horizontal motion and 
P is the wave period. This form was chosen to simulate the ef- 
fect of transverse motion of the flux sheet at the lower boundary. 
For simplicity we assume that all points of the lower boundary 
have this motion: this does not generate any waves in the ambi- 
ent medium, other than at the interface with the flux sheet. As 
a standard case in our simulation we use Vq = 750 m s"' and P 
= 24 s following Hasan et al. (2005). This short period is mo- 
tivated by the result of ( jHasan et al.,,2000 ) that high frequency 
waves would model the observational signature of wave heating 
less intermittent and thus in better agreement with the steady ob- 
served intensity from the magnetic network. We consider a uni- 
form horizontal displacement of the bottom boundary for half 
a period after which the motion is stopped (this corresponds 
to the impulsive case treated by| Hasan et al.|[2005| l. Such short 
duration motions are expected to be generated by the turbulent 
motion in the convectively unstable subsurface layers where the 



flux sheet is rooted. In terms of the analysis by Cranmer & van 



Ballegooijen (2005 1 of the kinematics of G-band bright points, 
this motion rather corresponds to a short, single step of their 
"random walk phase", for which these authors use a rms veloc- 
ity of 0.89 km s"' with a correlation time of bright-point motions 
of 60 s in accordance with the measurements of INisenson et al.l 
( 2003| l. The cases with higher velocities (see Table 2) would 
rather be representative of the "jump phase" for which |Cranmer 
& van Ballegooijen (2005 i use a velocity of 5 kms ' with a du 



ration of 20 s. This motion generates magnetoacoustic waves in 
the flux sheet. We first examine wave propagation and energy 
transport in a flux sheet with a sharp interface for the moderate 
and the strong field cases. In Sect. |6] we analyze the effect of 
varying the interface thickness. 



4. Dynamics 

4.1. Moderate field 

Let us consider a magnetic configuration in which the field 
strength at the axis of the flux sheet at z = is 800 G. In this 
case the /3 - I contour is well above the bottom boundary in 
the atmosphere and hence all the magnetic field lines emerging 
from the base of the sheet cross this layer at some height. Waves 
are excited at z = 0, where y8 > 1 (on the axis yS = 1.8), in the 
form of a fast (predominantly acoustic) wave and a slow (pre- 




400 800 
X(kni| 




Fig. 2. Temperature perturbations for the case in which the field 
strength at the axis at z = is 800 G (moderate field). The colors 
show the temperature perturbations at 40, 60, 80, and 120 s (from 
bottom to top) after initiation of an impulsive horizontal motion 
at the z = boundary of a duration of 12 s with an amplitude of 
750 ms"' and a period of P = 24 s. The thin black curves are 
field lines and the white curve represents the contour of (3 - 1 . 



dominantly magnetic|^ wave, which propagate respectively at 
the sound and the Alfven speeds. On the sheet axis, the acoustic 
and Alfven speeds at z = are 7.1 and 6.0 kms"', respectively 



' For brevity we call modes in the following simply acoustic and 
magnetic depending on the predominance of the thermal and magnetic 
nature of their restoring forces. 
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Fig. 3. Velocity components for the case in which the field strength at the axis at z=0 is 800 G (moderate field). The colors show the 
velocity components (a) V,,, along the field, and (b) V„, normal to the field, at 40, 60, and 80 s (from bottom to top) after initiation of 
an impulsive horizontal motion at the z = boundary of a duration of 12 s with an amplitude of 750 m s ' and a period of P = 24 s. 
The thin black curves are field lines and the white curve represents the contour of yS = 1. The field ahgned and normal components 
of velocity are not shown in the regions where B < 50 G. 



(see Table [T]l. The fast wave is created due to compression and 
rarefaction of the gas at the leading and trailing edge of the flux 
sheet, respectively: this can be clearly discerned in the snapshots 
of the temperature perturbation, AT (the temperature difference 
with respect to the initial value), shown in Fig. |2]at 40, 60, 80 
and 120 s after start of the perturbation. (These panels and pan- 
els in the following figures do not show the full height range of 
the computational domain but up to 1280 km only.) The black 
curves denote the magnetic field lines and the white curve de- 
picts the yS = 1 contour. The perturbations are 180° out of phase 



on opposite sides of the sheet axis. As these fast waves travel 
upwards they eventually cross the layer of yS - 1, where they 
change their label from "fast" to "slow", without changing their 
acoustic nature: this corresponds to a "mode transmission" in 
the sense of Cally (2007 1. The transmission coefficient depends 
(among others) on the "attack angle" i.e., the angle between the 
wave vector and the local direction of the magnetic field (Cally 
2007 ). On the /3 - 1 layer, away from the sheet axis, where the 
wave vector is not exactly parallel to the magnetic field, we do 
not have complete transmission of the fast wave to a slow wave. 
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Fig. 4. Temperature Perturbations for the case in which the field 
strength at the axis at z=0 is 1600 G (strong field) for times t - 
40, 60, 80, and 120 s. The coding corresponds to that of Fig.|2] 



Rather, there is a partial conversion of the mode from fast acous- 
tic to fast magnetic, so that the energy in the acoustic mode is 
reduced correspondingly. 

Figures[3^ and|3]3 shows the velocity components in the flow 
parallel (V,) and perpendicular {V„) to the field, respectively. The 
velocity components are shown only in regions where the field is 
greater than 50 G since in the ambient medium with weak field 
this decomposition is no longer meaningful. In general the waves 
possess both longitudinal and transverse velocity components, 
but in regions where y6 < 1, the parallel component essentially 



corresponds to the slow (acoustic) wave that is guided upward 
along the field. This correspondence can be seen by comparing 
the parallel flow pattern (in Fig. [3^) with the temperature pertur- 
bation in Fig.|2] 

The excitation at the bottom boundary also generates a slow 
(magnetic) wave with velocity perturbations normal to the mag- 
netic field line. In order to visualize the slow wave, we show 
the velocity component normal to the magnetic field in Fig.[3|3. 
The slow wave also encounters the layer of yS = 1 and undergoes 
mode transmission and conversion. Above the layer ofp- 1, the 
transmitted wave is a fast mode, which rapidly accelerates due 
to the sharp increase in Alfven speed with height. 



4.2. Strong field 

We now consider the case in which the field strength on the sheet 
axis is 1600 G (at z = 0). Here, the contour of /? = 1 approxi- 
mately traces the boundary of the flux sheet. The transverse mo- 
tion of the lower boundary generates slow (essentially acoustic) 
and fast (essentially magnetic) waves. Since the contour of /? = 1 
runs along the boundary of the flux sheet, waves generated in the 
sheet that travel upwards do not encounter this layer and hence 
do not undergo mode conversion. Figure |4] shows the tempera- 
ture perturbation AT at 40, 60, 80, and 120 sj^ Figure [s] shows 
the parallel and perpendicular components (with respect to the 
magnetic field) of the velocity. 

The slow (acoustic) wave is guided upwards along the field 
without changing character. On the other hand, the fast wave, 
which can travel across the field encounters the /3 - I con- 
tour at the boundary of the flux sheet. As the fast wave crosses 
this layer, it enters a region of negligible field and hence gets 
converted into a fast (acoustic) wave. This can be easily seen 
in the snapshot of temperature perturbations at an elapsed time 
of 40 s. The fast wave in the low-/? region, which is essentially 
a magnetic wave, undergoes mode conversion and becomes an 
acoustic wave, which creates fluctuations in temperature visible 
as wing like features in the periphery of the flux sheet between 
approximately z - 200 to 500 km. The fast wave gets refracted 
due to the gradients in Alfven speed higher up in the atmosphere. 
Furthermore, similar to Hasan et al.| ( |2005) , we find that the in- 
terface between the magnetic flux sheet and the ambient medium 
is a remarkable source of acoustic emission. It is visible in Fig.|4] 
as a wave of shell-like shape in the ambient medium that em- 
anates from the base of the flux sheet and subsequently propa- 
gates, as a fast acoustic wave, laterally away from it. 

Incidentally, the phase of transverse movement changes by 
180° between the moderate and strong field case as can be seen 
comparing Fig.jSJ) with Fig.[5]3. This is due to the development 
of a vortical flow from the high pressure leading edge of the 
flux sheet to the low pressure trailing edge that develops in the 
high-j6 photospheric layers of the moderately strong flux sheet 
but is largely suppressed in the strong field case, where it is from 
the beginning preceded by the fast (magnetic) wave that emerges 
right from the initial pulse. The development of a vortical flow in 
the moderate field case was also noticed in [Hasan et al.| ( |2005] l. 



- The temperature perturbations along the flux-sheet edges in the 
wake of the slow acoustic wave (red and blue ridges along the left and 
right boundary in the lower part of the flux sheet, respectively) do not 
pertain to a traveling wave. They are due to the finite shift of the flux 
sheet with respect to the initial, static configuration. This shift is com- 
pensated for by a corresponding shift of the unperturbed solution for the 
computation of energy fluxes in Sects. |5]and|6] 
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Fig. 5. Velocity components for the case in which the field strength at the axis at z-0 is 1600 G (strong field) for times t - 40, 60, 
and 80 s. The coding coiTesponds to that of Fig. [3] 



Besides the fast and slow acoustic and the fast magnetic 
wave that emanate directly from the initial perturbation there is 
also a slow magnetic mode from this source, which propagates 
in the high-y6 surface layer of the flux sheet. It is visible in Fig.jSj) 
as the yellow/red crescent-shaped perturbation, which trails the 
red and blue crescents pertaining to the fast magnetic mode. 
Different from the latter, which is maximal at the flux-sheet 
axis, the slow mode has maximal amplitude in the weak-field 
boundary-layer of the flux concentration. A similar but acoustic 
slow surface mode was found by Kho menko et al.| (j2008 ) when 
the driver was located in the high-/? layers of the flux concentra- 
tion. Here, this mode generates a remarkable amount of acoustic 
emission to the ambient medium as will be seen in Sect. |6l 



In a three-dimensional environment there would in general 
appear a third, intermediate wave type, additional to the slow and 
fast mode, depending on the geometry of the magnetic field and 
the initial perturbation. CoiTespondingly, we may expect mode 
coupling between all three wave types. In the presence of gravi- 
tational stratification, the y6 = 1 surface (more precisely the sur- 
face of equal Alfven and sound speed) would still constitute the 
critical layer for mode coupling so that we could expect scenar- 
ios not radically different from but more complex than those of 
Secs.l4l]andl4T2l 
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5. Energy transport 

We now consider the transport of energy in the various wave 
modes. Using the full nonlinear expression for the energy flux, 
it is not easily possible to identify the amount of energy caiTied 



by the magnetic and acoustic waves. Following Bogdan et al 
( |2003| l, we instead consider the wave flux using the expression 



given by Bray & Loughhead ( 1974 1 that represents the net trans- 
port of energy into the atmosphere: 



^PV + T-(^o ■ AB)V - -^(V ■ AB)Bo. 
An An 



(11) 



The first term on the right hand side of the equality sign is the 
net acoustic flux, and the last two terms are the net Poynting 
flux. The operator A gives the perturbations in the variable with 
respect to the initial equilibrium solution and Z?o refers to the 
unperturbed magnetic field. 

Figures [6^ and [S]? show the magnitude of the acoustic (left 
panels) and the Poynting flux (right panels) at 40, 60, and 80 s 
(from bottom to top) for the moderate field case. Since in the 
ambient medium the field strength is weak, the Poynting fluxes 
are not shown in this region. The Poynting flux is essentially 
the wave energy that is carried by the magnetic mode, which as 
expected, is localized to the flux sheet. On the other hand, the en- 
ergy transport in the acoustic-like component is more isotropic. 
At f = 40 s, we find from Fig.|6^ that the wave has just crossed 
the fi = I contour. Thereafter, it propagates as a slow wave 
guided along the field at the acoustic speed within the flux sheet 
and as a fast spherical-like wave in the surrounding quasi field- 
free medium. Inside the flux sheet, the energy in the magnetic 
component (Poynting flux) and the acoustic component is of the 
same order of magnitude. 

A comparison of these results with those for the strong field 
case (Fig. |7]i shows that in the latter case energy is transported 
by the fast wave much more rapidly, especially in the central 
regions of the flux sheet. This is due to the sharp increase of the 
Alfven speed with height above z > 200 km. At f = 40 s we 
find that the wave front associated with the magnetic component 
has already reached a height of about 500 km (close to the sheet 
axis), while the acoustic wave reaches this level only at about 
f = 80 s. 

From the contour plots of Figs.|6]and|7] we see that the fluxes 
in the ambient medium for the strong field case is still close to 
10'^ erg cm"^ s"', while for the moderate field, it is almost an or- 
der of magnitude less, suggesting that the flux sheets with strong 
fields are a more efficient source of acoustic fluxes into the am- 
bient medium. The "mode transmission" from fast (acoustic) to 
slow (acoustic) that takes place in the case of a moderate field, 
as explained in Sect. |4] can be seen in Fig|6^. Since the "at- 
tack angle" in this case is close to zero, a significant amount 
of acoustic transmission takes place across the layer of fi - 1. 
Another feature that we see in the plots of wave-energy fluxes is 
the "mode conversion" that takes place in the strong field case. 
The fast magnetic wave, which is generated inside the flux sheet 
can travel across the magnetic field. This mushroom like shape, 
which is expanding, can be easily discerned in the 40 s snapshot 
of the Poynting flux shown in Fig.[7j3. As this wave crosses the 
fi - I contour, it is converted into a fast (acoustic) wave. The 
wing like feature that can be seen in the 60 s snapshot of the 
acoustic fluxes (Fig. [7^) are due to the fast waves that have a 
moment ago undergone a "mode conversion" from magnetic to 
acoustic. 

Next we consider a field line to the left of the flux sheet axis, 
which encloses a fractional flux of 50%. The field aligned and 



the normal component of the wave-energy fluxes are calculated 
along this particular field line. Fig. [8] shows the positive, field 
aligned component of acoustic flux for the moderate and strong 
field case as a function of time and spatial coordinate z along the 
field line. The dotted curves in the figure show the space time po- 
sition of a hypothetical wavefront that travels with Alfven speed 
(steeper slope) and sound speed along this magnetic field line. 
With the help of this plots it is easy to separate the energy fluxes 
that is carried by the slow and the fast wave modes. The evo- 
lution of the fi - I layer is shown for the moderate field case. 
The perturbation of this layer as the wave crosses it can be seen 
clearly around 40 seconds. It moves down due to the decrease in 
pressure caused by the rarefaction front and then moves up when 
the compression front passes it. Most fraction of the flux lies 
parallel to the line that corresponds to the hypothetical acoustic 
wave, which is a slow mode in the region where yS < 1 . In the 
strong field case, above approximately z - 800 km and for times 
t S 120 s, the acoustic flux carried in the compressive (trailing) 
phase starts to catch up the slightly slower moving expansive 
phase and the flux gets confined into a narrow shock forming 
region. This is also visible in the case of the moderately strong 
field. This behavior is not present along the coiTesponding field 
line to the right of the sheet axis (not shown here), where the 
compressive phase is leading so that the compressive and expan- 
sive phase of the perturbation slightly diverge with time. 

The acoustic fluxes are of the order 10^ ergs cm"^ s"^ The 
Poynting fluxes caiTied by the fast mode in this region can be 
identified by the coloured contours that gather along the dotted 
curves corresponding to the hypothetical Alfven wave (Fig. |9]l. 
Comparing the two fluxes (Fig. [9] with Fig.|8]l, it is clear that the 
acoustic flux caiTied by the slow mode is larger than the Poynting 
flux, especially in the moderate field case. The Poynting flux 
rapidly weakens with time and height because it is not guided 
along the field lines like the slow mode but rapidly diverges 
across the field and part of the Poynting flux converts to acoustic 
again as explained in Sect.|4] Also from Fig.|9]it can be seen that 
while the magnetically dominated fast mode starts right from 
the excitation level at z = in the strong field case, it starts in 
the weak field case only after about 40 s when the fast acoustic 
wave reaches the conversion layer where fi ^ \ and partially un- 
dergoes mode conversion. Therefore, the fast (magnetic) mode 
is from the beginning weaker in the moderate as compared to 
the strong field case. 

Table |2] shows the temporal maximum of the horizontally 
averaged vertical components of acoustic and Poynting fluxes at 
three different heights for the strong field case. We have con- 
sidered three different amplitudes and periods for the initial ex- 
citations. Although the field aligned acoustic fluxes on the spe- 
cific field line considered in Fig. [8] reach values of the order of 
10^ erg cm"^ s"' at a height of z=1000 km, the horizontally av- 
eraged fluxes are typically an order of magnitude less, depend- 
ing upon the amplitude of the initial excitation. The Poynting 
fluxes shown in the table are the maximum value that the fluxes 
reach in the interval between the start of the simulation until the 
time when the fast wave reaches the top boundary (around 60s). 
Hence these fluxes are due to the fast mode for z - 500 km and 
z - 1000 km, since within this time limit the slow mode has not 
yet reached these heights. The Poynting fluxes associated with 
the fast mode are relatively lower in magnitude compared to the 
acoustic fluxes. It should be noted that there is also considerable 
Poynting flux associated with the slow mode, since these waves 
also perturb the magnetic field. The acoustic fluxes of the mod- 
erate field case reach only less than 70% of that of the strong 
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X (km) X (km) 

a) Acoustic b) Poynting 

Fig. 6. Wave-energy fluxes (absolute values) for the case in which the field strength at the axis at z = is 800 G (moderate field). The 
colors show (a) the acoustic flux, and (b) the Poynting flux, at 40, 60, and 80 s (from bottom to top) after initiation of an impulsive 
horizontal motion at the z = boundary of a duration of 12 s with an amplitude of 750 m s"' and a period of f = 24 s. The thin 
black curves are field lines and the thick black curve represents the contour of /3 - 1 . The Poynting fluxes are not shown in the 
ambient medium where B < 50 G. 

Table 2. Temporal maximum of the horizontally averaged, vertical component of the wave-energy fluxes for the strong field case 



Initial Excitation 


Fa., 


(10^ erg cm" 






(10* erg cm" 


s-') 


z=100km 


z=500 km 


z= 1000 km 


z=100km 


z=500 km 


z= 1000 km 


0.75 km s"', 


24s 


11.36 


1.96 


1.33 


29.38 


1.08 


0.14 


0.75 km s-', 


120s 


35.75 


27.70 


4.02 


134.29 


0.79 


0.07 


0.75 km s-', 


240s 


20.90 


8.58 


3.30 


131.84 


0.36 


0.02 


1.50 kms"', 


24s 


44.55 


7.68 


3.34 


115.79 


4.29 


0.57 


3.00 km s"', 


24s 


16841 


30.40 


6.22 


434.03 


16.90 


2.31 
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Fig. 7. Wave-energy fluxes for the case in which the field strength at the axis at z = is 1600 G (strong field) for times t - 40, 60, 
and 80 s. The coding corresponds to that of Fig.|6] The Poynting fluxes are not shown in the ambient medium where B < 200 G. 



field configuration and the Poynting fluxes are negligible in in 
this case. 



6. Effects of the boundary-layer width 

We now study the acoustic emission of the magnetic flux con- 
centrations into the ambient medium by varying the width of 
the boundary layer between the flux sheet and ambient medium. 
This is caiTied out by comparing the result of simulations with a 
sharp interface of width 20 km to that with a width of 80 km (see 
Fig. [T]i, where the width can be varied by choosing appropriate 
values of and 1^2 in Eq. (|2]i. 



We examine the acoustic emission from the two peripheral 
(control) field lines to the left and to the right of the flux-sheet 
axis that encompass 90% of the magnetic flux. These correspond 
to the outermost field lines that are plotted in Figs.|2]to|7] These 
field lines are located in the high-/? region with fi » 1, all the 
way from the base to the merging height, where the flux sheet 
starts to fill the entire width of the computational domain. The 
acoustic emission from the peripheral field line to the right and 
to the left of the flux-sheet axis is practically identical. 

Figure [TO] shows the acoustic emission from the flux sheet 
into the ambient medium for the peripheral field line to the left 
of the flux sheet with the strong field (Bq - 1600 G) and the 
cases of the sharp interface (left panel) and the wide interface 
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Fig. 8. The field aligned positive (upwardly directed) component of acoustic wave-energy flux as a function of time on a field line 
on the left side of the axis that encloses a fractional flux of 50%. The left panel represents the case in which the field strength at the 
axis at z = is 800 G (moderate field) and the right panel represents the strong field case (1600 G). 
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Fig. 9. The field aligned positive (upwardly directed) component of the Poynting flux as a function of time on a field line on the left 
side of the axis that encloses a fractional flux of 50%. The left panel represents the case in which the field strength at the axis at 
z = is 800 G (moderate field) and the right panel represents the strong field case (1600 G). 



(right panel). Concentrating on the case with the sharp interface 
first, we see that acoustic flux is initially generated by the fast 
mode that stems from the transversal motion of the flux sheet to 
the right hand side, which causes a compression and expansion 
to the right and left side of the flux-sheet edge, respectively. This 
movement generates a net acoustic flux away from the flux sheet 
on both sides. It is visible in Figs. |4]and|7]as the shell-like an- 
tisymmetric wave that emanates from the base of the flux sheet 
propagating into the ambient medium. At a height of z = 100 km 
the peak value of this flux is 3 x 10** erg cm"^ s"' for the sharp in- 
terface but only 1 .2x 10^ erg cm"^ s"' for the wide interface. This 
is because the sharp interface acts like a hard wall that pushes 
against the ambient medium, while the wide interface is more 
compressible and acts more softly. 

Near the flux-sheet boundary this wave seamlessly connects 
to the tips of the crescent-like fast (magnetic) mode of the flux- 
sheet interior as can be best seen when comparing the first two 
snapshots of Figs. |4] and |5]5. There, acoustic flux is generated 
by continuous leakage and conversion from the magnetic mode, 
giving rise to the steeper of the two horizontally running, in- 
clined ridges of acoustic flux, visible in the lower part of both 
panels of Fig. [TO] This leakage is more efficient in the case of 
the wide interface than in the case to the sharp interface so that 
the corresponding ridge extends over a longer period of time in 
the former compared to the latter case. However, it cannot com- 



pensate for the larger initial flux that emanates from the more 
confined (sharp) boundary. 

Starting at about f = 25 s in case of the sharp interface, one 
can see a less steep and weaker branch of acoustic flux that is 
connected to the slow (magnetic) mode that propagates in the 
high-j6 boundary-layer of the flux sheet. Obviously it creates a 
non-negligible source of acoustic flux to the ambient medium. It 
is also present in case of the wide interface. 

The two horizontally running ridges of acoustic flux in the 
case of the sharp interface (Fig. 



10 



left) is slightly more in 
clined compared to the case with the wide interface (Fig. 10 
left), where the peripheral (control) field line expands more in 
the horizontal direction so that the wave travels a longer distance 
to reach it. 

At approximately f = 45 s we start to see acoustic flux ap- 
pearing at a height of about z = 1000 km. This flux originates 
from the refracted fast (magnetic) wave within the flux sheet. 
Since this wave quickly accelerates and refracts with height, it 
reaches the flux-sheet boundary sooner at z = 1000 km than in 
the height range 500 km < z < 800 km. This wave undergoes 
conversion from fast, predominantly magnetic to fast, predom- 
inantly acoustic as it crosses the region where yS = 1. Because 
it travels essentially perpendicular to the field near the flux- 
sheet boundary, the conversion is particularly efficient. While 
this ridge of acoustic flux originates from the leading phase of 
the fast (magnetic) wave that corresponds to a movement to the 
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Fig. 10. Acoustic flux perpendicular to the peripheral field lines that encompass 90% of the magnetic flux as a function of time 
and height along the field line. Only the outwardly directed flux is shown. Left: strong field case with a sharp interface between 
flux-sheet interior and ambient medium. Right: strong field case with a wide interface. 

Table 3. Total acoustic emission from the flux sheet into the ambient medium for different boundary layer widths. 



Boundary-layer width 


Total acoustic emission 




(10'* ergcm"') 


20 km (sharp) 


23.40 


40 km (medium) 


13.22 


80 km (wide) 


8.13 



right (red big crescent in the 40 s snapshot of Fig. [S]?), a second, 
parallel running negative ridge, stems from the following phase, 
corresponding to a movement to the left (blue crescent in the 40 s 
snapshot of Fig. 

Table |3] shows the total acoustic emission to the ambient 
medium, still from and perpendicularly across the field lines that 
encompasses 90% of the total magnetic flux for cases with 3 
different boundary layer widths. The energy is computed by in- 
tegration of the perpendicular flux along the peripheral control 
field lines to the left and to the right over the full height range of 
the computational domain and from f = s to f = 64 s for unit 
width. The total acoustic energy leaving the flux sheet with the 
wide interface is only 35% of that with the sharp interface. In 
this sense, a flux sheet with a sharp interface is more efficient in 
providing acoustic flux to the ambient medium than a flux sheet 
with a wide interface as conjectured by Hasan et al. (2005 i 



7. Summary 



This work is an extension of the previous work done by |Hasan| 
et al. ( |2005[ l. Wave excitation occurs in a magnetic flux con- 
centration by a transverse motion of the base. The present work 
extends the previous calculations to the case of a flux sheet with 
moderate field strength. In addition, a new feature of the present 
work is that we estimate the energy carried by the waves and we 
examine the effect of varying the thickness of the tube-ambient 
medium interface on the acoustic emission. 

We have found that the nature of the modes excited depends 
upon the value of /3 in the region where the driving motion oc- 
curs. When j6 is large, the slow wave is a transverse magnetic 
mode that propagates along the field lines and undergoes mode 
transmission as it crosses the /? = 1 layer. In this case, the wave 
only changes label from slow to fast, but remains magnetic in 
character throughout the flux sheet. The fast mode, which prop- 
agates almost isotropically, undergoes both mode conversion and 



transmission at the y6 = 1 surface depending on the "attack an- 
gle", the angle between the wave vector and the magnetic field. 
On the other hand, in the case of a strong magnetic field (low-/? 
case), where the level of /3 - 1 is below the driving region, the 
fast (magnetic) and slow (acoustic) modes propagate through the 
flux sheet atmosphere without changing character. 

We find that the magnetically dominated fast wave within 
the low-yS region of the flux sheet undergoes strong refraction so 
that it finally leaves the flux sheet in the lateral direction, where it 
gets partially and mainly converted to a fast, acoustically domi- 
nated wave. This eff'ect is particularly visible in the case of a flux 
sheet with strong magnetic field. 

We also see an asymmetry in the wave structure on both sides 
of the flux sheet axis. This comes because the leading front of 
the predominantly acoustic mode is compressional on the one 
hand side and expansive on the other side and vice versa for the 
following phase. Since the compressive phase travels faster as 
the sound speed is larger, the two phases move either apart from 
each other or converge. In principle, this asymmetry should give 
rise to observable signatures. 

Recent observations of the chromospheric network are sug- 
gestive of Can network grains associated with plasma with 
quasi-steady heating at heights between 0.5 and 1 Mm inside 
magnetic flux concentrations ( [Hasan & van Ball egooijen 2008J. 
Let us now estimate the acoustic energy flux transported into 
the chromosphere through a single short duration pulse as has 
been treated in the present paper (the magnetic modes are almost 
incompressible and not efficient for heating the atmosphere). 
We consider the energy flux at a height of 1000 km. For the 
strong and moderate field cases, the maximum values of the 
acoustic fluxes at z = 1000 km are ~15xlO^ erg cm"^ s"' and 
~5xlO^ erg cm"^ s"', respectively. However, it should be noted 
that although the fluxes can reach values up to 10^ erg cm"^ s"', 
the spatially averaged values are much less. From Table |2] we 
obtain for the strong field case a temporal maximum of the hor- 
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izontally averaged acoustic flux at z = 1000 km of a few times 
10^ erg cm"^ s depending on the excitation amplitude and pe- 
riod. 

We see that the strong field configuration is a more efficient 
source of acoustic waves in the ambient medium compared to the 
weak to moderate field configurations. For the cases considered 
here, they differ by almost a factor of two. The width of the tran- 
sition layer between the flux sheet and the ambient medium has 
significant eff'ect on the acoustic wave emission as was initially 
conjectured by Hasan et al. (2005 1. Our quantitative calculations 
substantiate this hypothesis: a flux sheet with a sharp interface 
emits almost four times the energy emitted by a flux sheet with 
a wide interface. 



8. Discussion and conclusion 

The energy losses in the magnetic network at chromospheric 
heights are of the order of 10^ erg cm"^ s"'. Even though the 
acoustic energy flux produced by the transverse excitation move- 
ment can temporarily reach this value at certain locations, the 
values of Table |2] tell us that in the spatial average these en- 
ergy losses cannot be balanced by the acoustic energy flux gen- 
erated in our model. This conclusion is emphasized by the fact 
that the values of Table |2] are temporal maxima: the temporal 
mean would be lower. In order to be compatible with the ob- 
served quasi-steady Ca emission the injection of energy needs 
to be in the form of sustained short duration pulses as argued by 
( Hasan & van Ballegooijen|2 008 ) but these pulses could proba- 
bly not maintain the maximum values of acoustic flux as quoted 
in Table 12] 

Possibly with the exception of the case corresponding to the 
last row in Table |2] the transverse excitation considered here 
rather correspond to the "random walk phase" of the model by 
Cranmer & van Ballegooijen (2005]l. Excitations corresponding 
to the "jump phase" with even higher velocity amplitudes than 
considered in the present paper might temporarily be capable of 
providing the required energy flux. However, with a mean inter- 
val time of 360 s these jump events are probably not responsible 
for the heating observed in Ca ii network grains, which requires 
a more steady or high frequency source. 

We have not considered photospheric radiative losses, which 
would considerably damp the waves before they reach chromo- 
spheric heights (Carlsson & Stein 2002). If these radiative losses 
are taken into account, the fluxes would further be lowered. Also 
not all of the acoustic energy flux would be available for radia- 
tive energy loss in the chromosphere depending on the details of 
this NLTE process. All this implies that acoustic waves gener- 
ated by transverse motions of the footpoints of magnetic network 
elements cannot balance the chromospheric energy requirements 
of network regions. 

This conclusion cannot be expected to drastically change 
when turning to three spatial dimensions. The details of the 
mode coupling and the partition of energy fluxes to the various 
modes would become more complex but the share of energy that 
resides in the acoustic mode cannot be much larger than in the 
two-dimensional case. On the contrary, the energy flux gener- 
ated at the footpoint of the magnetic element would have to be 
distributed to a larger area in three spatial dimension so that the 
spatial mean at z = 1000 km would be lower. 

We have only considered single, short duration, transverse 
pulses for the wave excitation. A more realistic driver with sus- 
tained pulses of varying lengths, velocities, and time intervals 
would give rise to highly non-linear dynamics, which might 
yield increased acoustic fluxes. Also we have not considered 



longitudinal wave excitation, which would be available primar- 
ily from global p-mode oscillations. The latter are expected to 
provide low frequency slow mode waves to the outer atmo- 
sphere via magnetic portals in the presence of inclined strong 
magnetic fields, where they would be available for dissipation 
through shock formation ( Michalitsanos| 1 973 Suematsu|1990 



[Hansteen et al.|2006) [Jefferies et al.|2006| l. In fact, this mecha 
nism would also work in the periphery of a vertically oriented 
flux tube, where the field is strongly inclined with respect to the 
vertical direction. Another source of energy that was not consid- 
ered here may come from direct dissipation of magnetic fields 
through Ohmic dissipation. 
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